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Abstract 

The equation for the charge vertex 7 of the t — J model is derived and solved 
in leading order of an 1 /N expansion, working directly in terms of Hubbard 
operators. Various quantities which depend crucially on 7 are then calculated, 
such as the life time and the transport life time of electrons due to a charge 
coupling to other degrees of freedom and the charge-charge correlation func- 
tion. Our results show that the static screening of charges and the dynamics 
of charge fluctuations depend only weakly on J and are mainly determined 
by the constraint of having no double occupancies of sites. 
PACS numbers: 72.12.Di, 63.20.Kr., 74.72.-h 
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The charge fluctuation spectrum and the screening properties of the t-J model are largely 
determined by the charge vertex 7(fc,g). Here, k denotes both the momentum k of an 
electron and the Matsubara frequency iun- Similarly, q stands for the transferred momentum 
q and the frequency A simple interpretation of 7 has been given in0: It multiplies the 
bare electron-phonon interaction g{k,q) to yield an effective electron-phonon interaction 
which takes into account screening effects due to the constraint of no double occupancies of 
sites. Another property of 7 is that its poles in the second frequency argument determine 
the dispersion of collective density waves and thus is an important ingredient for the density- 
density correlation function D{q). For the case J = the properties of the static vertex 
7 has been investigated ini'i. In particular, it has been shown that 7 exhibits a strong 
momentum dependence in q at low frequencies and small and intermediate dopings which 
tends to suppress the effective electron-phonon interaction in the tt'-model {t amd t' denote 
hopping integrals between nearest and second- nearest neighbor sites, respectively). This 
effect is especially pronounced in transport quantities. The frequency dependence of 7 and 
D has been investigated, again for J = 0, ini. There it has been shown that D is nearly 
exclusively determined by collective effects and has an energy scale substantially larger than 
the effective band width, in agreement with computer simulations^. The purpose of this 
communication is to extend the above results to the t-J model and to investigate to what 
extent the above properties of 7 and D depend on J. 

The Hamiltonian of the t-J model reads 

H = - %xfxf+ jk^rxf. (1) 

p— l...A^ p,q— i...Ar 

The subscripts i, j stand for lattice sites; the superscripts p, q denote for p = the unoccupied 
and for p = 1...N singly occupied states with a spin index p. This means that the original 
SU (2) spin space has been extended to a SU{N) space which is a well-known procedure in 
slave boson calculationJl. The Hubbard operators Xf^ with p = 0, g = and p > 0, q > 
have bosonic and those with p = 0,g>Oorp>0,g = have fermionic character. They 
obey the following commutation and anticommutaion rules, respectively. 



[Xf\ ]^ = 5,,{5,rXr T Ssr,Xn- (2) 

In the SU{N) model, considered here, the X-operators are subject to the constraint 

A N 

E^r = y (3) 

This means that at most N/2 of the N states at each site can be occupied at the same 
time. The first term in Eq.(l) describes the hopping of particles between the sites i and 
j with matrix elements tij. The second term in Eq.(l) denotes the Heisenberg interaction 
between the spin densities at site i and j with the exchange constants Jij. In the following 
we consider Jij only between nearest neighbors {Jij = J) and tij between nearest {tij = t) 
and next nearest {tij = t') neighbors. The coupling constants in Eq.(l) have been scaled 
with N in such a way that the limit N ^ oo describes an interesting physical case and that 
for N = 2 the usual t-J model is, except for an overall factor 1/2, recovered. 

Using a 1/N expansion the Hamiltonian Eq.(l) has been investigted for J = in0 
and, in more detail, ini. These treatments can be generalized to a finite value of J in a 
straightforward way: The equation of motion for a fermionic X-operator is, using Eq.(l): 

^X°«(l)= J2 /f^2c/3t(°^'^'^'°^')X^'2«2(2)X°^«(3) (4) 
^'''i P2,g2,g3"' 12 3 

with 

^f) = 5(1 - 2)6,,oKo^.,Al - 3)/X 

+ ^.,sK,.m - mi - 3) - 6{1 - 3) J(l - 2)/2)/X. (5) 

Here, 1 is an abbreviation for iiTi, i.e., 1 = (?iTi), where ri denotes the imaginary time. 
t{l — 2) is equal to ti-^i26{ri — T2) and J(l — 2) equal to Jij^i2S{Ti — T2). The first term in 
the parantheses on the right-hand side of Eq. (5) describes hopping without flip of the spin 
whereas the second one hopping with a spin-flip. Comparing the above Eqs.(4) and (5) with 
Eqs.(9) and (10) of Ref.i one finds that the Heisenberg term in H just adds a contribution 
to the spin-flip hopping term. The perturbation expansion ini rests on two relations: The 
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equation of motion for fermionic Hubbard operators and Eq.(31) ina which relates expec- 
tation values of bosonlike Hubbard operators to Green's functions. The first relation is 
modified by the Heisenberg term in the above way, the second relation is unchanged. As a 
result, it is straightforward to generalize the expressions for the self-energy, the vertex etc. 
ini to the case of a finite J. 

Using the above procedure one obtains from Eq.(37) ini the following expression for the 
self-energy in 0(1) of the t-J model: 

S(l - 2) = 6(1 -2) J d3t{l - 3)g{3+ - 1) - t{l - 2) < X''{1) > - 2). (6) 

The normalized Green's function g (denoted by G ini) is related to S via Dyson's equation 

I d3{ - b{\ - 3)— - E(l - 3))gi3 - 2) = 6(1 - 2). (7) 

< X^^{1) > is the expectation value of X°°(l). Both, E and g are diagonal in the in- 
ternal indices so we have omitted them in the above equations. The self-energy in Eq.(6) 
is instantaneous giving rise to a frequency-independent, but momentum-dependent renor- 
malized one-particle energy e(k). After a Fourier transformation Eqs.(6) and (7) yield 
g{k, iuri) = 1/ {iuJn - ^(k)) with 

e(k) = A - got(k) - ^ E -^(k + P)/(e(p)). (8) 

Here we have ^(k) = e(k) — /i and A = J2pt{p) f {(,{p)) , where / is the Fermi function, 
Nc the number of primitive cells and go = '^/2 with the doping 6. 

Taking the Heisenberg interaction also into account the vertex equation (39) ini becomes 

f(12; 3) = 5(1 - 2)5(1 - 3) + t(l -2) J dhdQg{l - 5)f(56; 3)^(6 - 1+) 
+5(1 -2) j dAdhdQt{l - A)g{A - 5)f (56; 3)^(6 - 1) - '^'^^ ~ j d5d6g{l - 5)f (56; 3)^(6 - 2). (9) 

Writing r(12; 3) = 7(1 — 2, 1 — 3) Eq.(9) becomes after a Fourier transformation 

7(A;, q) = l + jT- + t{k' + g) - ^ _ ' )9{k')g{k' + qh{k', q). (10) 



Since the k-dependence of J and t are given by trigonometric or products of trigonometric 
functions Eq.(lO) represents an integral equation with a kernel consisting of 6 separable 
contributions. Eq.(lO) thus can be reduced to a 6x6 matrix equation with the solution 

7(fc, g) = 1 - E E(^"/3 + x{q)Uxp2{q). (H) 

o=l 13=1 

XaM=T.Ga{k',q)F^{k'). (12) 

k' 

The vectors F and G are given by 

Fa(k) = (t(k), 1, Jcoskx, Jsinkx, Jcosky, Jsinky) (13) 

and 

Ga{k, q) = (1, t(k + q), cosk^, sinkx, cosky, sinky)Il{k, q) (14) 

with n(A;, q) = —g{k)g{k + q). The frequency sum in Eq.(12) involves only 11 and can easily 
be carried out: 

^n(.,,,^iMi±MM. (15) 

V ^(k ) - Uq + k') - ^^^n 

Calculation of 'y{k, q) thus requires essentially the calculation of the susceptibility matrix 
Xajsiq)- For J = the matrix inversion in Eq.(ll) can be done explicitly and one obtains 
Eqs.(12)-(15) oil. In contrast to this special case, 7(/c,g) depends for J 7^ also on the 
vector k for a given doping. 

The Green's function D{q) for density fluctuations is given by 

^(?) = ^E7(A:,g)n(A;,g). (16) 

k 



Using Eq.(ll) we obtain 



D{q) = -j2{l + x{q))TjxM. (17) 

13=1 



Carrying out the analytic continuation u + irj the density-density correlation function 

D{c{,uj) is equal to the negative imaginary part of D{q,uj + irf). 



We have evaluated numerically the susceptibility matrix Xapiq) using a typical mesh 
of 1000x1000 points in the Brillouin zone for the summation over k'. Fig.l shows results 
for the zero-frequency vertex 7(k, q) for J = 0.2, t' = (all energies are in units of the 
nearest-neighbor hopping energy t) and three different dopings 6. For each doping k is put 
on the Fermi line in (1, l)-direction and q is varied along the (1, l)-direction from zero to 
the maximal momentum transfer for points on the Fermi line. Fig.l should be compared 
with Fig.l of^ where a similar plot for 7 is given for the case J = t' = 0. (The momentum 
scale in that reference should be scaled by a factor in order to have the same absolute 
scale) . The Figures clearly show that the main property of the vertex found in ^ for J = is 
also present for J = 0.2: For large dopings 7 varies only smoothly with momentum whereas 
at smaller dopings 7 develops a forward scattering peak with a width ~ 6 due to a strong 
suppression of large momentum transfers. This implies that the effective charge interaction 
of electrons with other degrees of freedom (impurities, phonons etc.) is essentially the bare 
one at small but heavily suppressed at large momenta. The curve for 6 = 0.10 in Fig.l 
shows a new feature: It passes through zero at a small momentum which means that the 
effective interaction is exactly zero at this point due to correlation effects. With increasing 
momentum it goes through a minimum with a negative value and approaches zero from 
below at large momenta. A similar, but less pronounced behavior, has been foundi in the 
one-dimensional t-model and in the two-dimensional tt'-model with a finite t'. Finally, we 
have chosen in Fig.l a rather small value for J and not too small values for 6 in order to 
avoid singularities in 7 due to instabilities of the homogenous phase^. 

The above vertex function allows to answer the following question: How much are the 
inverse life time 1/r and inverse transport life time 1 / Ttr of an electron affected by electronic 
correlations if the coupling of the electron to additional degrees of freedom (phonons, impu- 
rities etc.) is due to the interaction of charges? The answer becomes especially simple if one 
assumes that the bare coupling function is structureless, i.e., is independent of momentum 
and frequency. Ini is has been shown that the quantities Ai,Atr defined by 
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C « 



7(k,k-k' 




>k>k', 



(18) 



1 — 



A. 



C « 



7(k,k-k' 




(v(k) - v(k'))' >k>k' /(2 « v2(k) >k>k', 



(19) 



describe changes in the inverse hfe time 1/r (or, in the Ehashberg function a'^F{uj) for s- 
wave supercondcutivity) and the inverse transport hfe time 1 / r^^ due to correlation effects. 
The overaU factor C is chosen such that Ai = Atr = 1 for 6-^1, i.e., the empty band hmit. 
If 7 in Eqs.(18) and (19) depends only weakly on momentum we have Ai ~ Atr- On the 
other hand, if 7 is nonzero only for k ~ k' Atr is much smaller than Ai. Fig. 2 shows Ai and 
Atr as a function of 6 for J = 0.3 and t' = —0.25. With decreasing doping Ai and Atr first 
pass through a maximum at around 5 ~ 0.8 and then decrease monotonically by around a 
factor 2 and 4, respectivley, until 6 ~ 0.2. The more and more pronounced appearance of 
a forward scattering peak in 7 at still smaller dopings would cause a further decrease in Ai 
and Afr and especially in the ratio Atr/ Ai. However, we exclude this low-doping region from 
our considerations because of the occurrence of instabilities of the homogenous phase in that 
region. Figure 2 suggests that correlation effects suppress Ai and, even stronger, Atr, moving 
from the overdoped towards the maximal doped regime. Fig. 10 ini presents results for Ai 
and Atr and for (using our energy units) J = and t' = —0.20. Comparing this Figure with 
our present Figure 2 one concludes that Ai and A^.^ depend only very weakly on J. (The 
additional interpretation of the quantity 6 ■ Atr as being proportional to the resistivity ini 
should be dropped since the Drude weight entering the static part of the resisitivity depends 
strongly on 6 and the resisitivity near half-filling is characterized by Atr/ 5 rather than by 



Fig. 3 shows the density-density correlation function -D(q, uj) for J = 0.1 (left panel) and 
J = 0.3 (right panel) for various momenta q. Curves corresponding to the same momentum 
practically coincide with each other demonstrating that D(q, u;) and thus also the dynamic 
part of the vertex are nearly independent of J. This implies, that D is dominated by 
collective effects in form of an infinitely sharp, dispersive sound peak also in the presence 



Atr-S). 
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of the Heisenberg interaction. This peak has been broadened in Fig. 3 by using a finite 
value of 0.1 for rj. The energy of this peak is in general much larger than the width of the 
renormalized band ( which is 0.96 for J = 0.1 and 1.28 for J = 0.3). The contribution 
of the particle-hole contimuum to D is nearly invisible if the sound peak is well above the 
particle- hole spectrum like, for instance, in the case q = (7r,7r). In the other cases like 
q = (27r/5, 7r/5), (tt/S, 37r/5), or(7r, 0) the collective peak is not so well separated from the 
particle-hole continuum and D has structure also at low frequency refiecting density of states 
of single particle-hole excitations. The absence of an noticeable dependence of the peak 
position on J as well as the quite different energy scales for charge and spin fluctuationJl 
remind of spin-charge separation found in one-dimensional models. 

In conclusion, the equation for the charge vertex 7 of the t — J-model has been derived in 
leading order of an 1 /N expansion, reduced to a 6x6 system of linear equations, and solved 
numerically. Our results for the momentum and frequency dependence of 7 show only a 
weak dependence on J. We also discussed various properties which depend sensitively on 7, 
namely, the effect of correlations on the inverse life time and the inverse transport life time 
of electrons and the dynamics of charge fluctuations. Our conclusion is that these quantities 
depend only weakly on J and are mainly determined by the constraint of having no double 
occupancies of sites. These flndings are consistent with recent exact numerical results from 
small cluster^'!. 

Acknowledgement: One of us (M.L.K.) would like to thank Prof. Michael Mehring 
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FIGURES 

FIG. 1. Zero- frequency vertex function 7(k, q) as a function of aq with k fixed on the Fermi 
Une along the (1, l)-direction for three dopings S. 

FIG. 2. Correlation induced enhancements Ai and A^^ as a function of doping S for J = 0.3 
and t' = -0.25. 

FIG. 3. Density-density correlation function D(q,Lo) as a function of energy lo for 6 = 0.2 and 
J = 0.1 (left panel) and J = 0.3 (right panel) using rj = 0.1. 
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